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Abstract 

Evolution of the main signal in a Lorentz dispersive medium is 
considered. The signal propagating in the medium is excited by a 
sine-modulated pulse signal, with its envelope described by a hyper- 
bolic tangent function. Both uniform and non-uniform asymptotic 
representations for the signal are found. It is shown when the uniform 
representation can be reduced to the non-uniform one. The results 
obtained are illustrated with a numerical example. 

1 Introduction 

Investigations on propagation of pulse signals in dispersive media date back 
to the beginning of 20th century. The fundamental research in this area is 
due to Sommerfeld [1] and Brillouin [2, 3]. Although steady interest in this 
kind of propagation was observed in the literature since then, a new impetus 
has been added recently due to new applications of the theory in fiber-optics 
communication and integrated-optics. Also, the knowledge of pulse propa- 
gation in a dispersive medium, and of accompanying electromagnetic energy 
losses in the medium, became of vital importance in radiotherapy. A signifi- 
cant contribution to the research on dispersion phenomena in Lorentz media 
is due to Oughstun and Sherman [4]. Equipped with better computation 
techniques and advanced asymptotic methods, they extended the analysis 
to models more closely reflecting practical applications. In particular, they 
considered signals with a finite rise time and employed uniform asymptotic 
expansions in their analysis. (Uniform expansions remain valid as their pa- 
rameters vary while non-uniform expansions break down at some parameters 
values.) 

Here, we also consider the evolution of the main signal excited in a Lorentz 
dispersive medium by a signal with a finite rise time. However, unlike Ough- 
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stun and Sherman work, where the envelope of the initial signal is described 
by an everywhere smooth function of time which tends to zero as time goes 
to minus infinity, our exciting signal is switched abruptly at a finite time 
instant, and vanishes identically for earlier times. In the analysis carried out 
in this paper we apply the Bleistein and Handclsman [10] theory of uniform 
asymptotic evaluation of integrals with nearby saddle point and an algebraic 
singularity. We show, how the uniform representation of the evolution of the 
main signal reduces to the non-uniform representation, which can otherwise 
be obtained by residues. 

The results obtained here are illustrated with a numerical example. 



2 Formulation of the problem 

Consider the problem of an electromagnetic plane wave propagation in a 
homogeneous, isotropic medium, whose dispersive properties are described by 
the Lorentz model of resonance polarization. The complex index of refraction 
in the medium is given by [4] 

Here, 6^ = A7rNe^/m, where N, e and m represent the number of electrons 
per unit volume, electron charge and its mass, respectively, 5 is a damping 
constant and ujq is a characteristic frequency. 

Let the signal /lo(t) in the plane 2; = be a sine wave of a fixed real 
frequency cUc with its envelope described by a real function u{t), identically 
vanishing for t < 0, i.e. 

. . . _ / t<0 , . 

u{t) smM t>0. 

Then arbitrary component of the wave propagating in the direction of in- 
creasing z (or of a corresponding Hertz vector) can be represented in the 
medium by the scalar function [4] 

A{z,t) = — Re < i / M(a; — a;c) exp -(f)(uj,6) duj>, (3) 
2n I Jia-00 ic J J 

where u{u!) is the Laplace transform of u{t). The complex phase function 
(f){uj,6) is given by 

(t){uj,e)^iuj[n{uj)-ei (4) 
where the dimensionless parameter 

.^f (5, 
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describes the space-time point {z, t). 

It is here assumed that the envelope of the incident pulse is described by 

= { t°anh/5t l>l 

where the parameter /3 > determines the rate of the pulse growth. 
The Laplace transform of u{t) is 

and the beta function B is related to the psi function t/j by [9] 




By using (7) in (3) we obtain the formula 

A{z, t) = —Rc\i Ufs{u; - uj,)e-c'>'^'^''^^ du \ , (9) 

which describes the dynamics of the signal excited at 2; = by Ao{t), and 
propagating in the Lorentz dispersive medium in the direction of growing z. 
The uniqueness of this solution is proved in Sec. 2 of [3] . 

In this work we study the poles contribution to the asymptotic expansion 
of A{z, t). We denote this contribution by Ac{z, t) and find both non-uniform 
and uniform asymptotic expressions for it. 



3 Non-uniform asymptotic expression for 

Ac{z, t) 

In finding an asymptotic expansion of the integral defined by (9) it is essen- 
tial to determine the location of its critical points, including saddle points 
and the poles in the complex a;-plane. The equation governing the location 
of the saddle points does not seem to be solvable exactly. Instead, different 
approximate solutions were obtained by Brillouin [3], Kelbert and Sazonov 
[6], and Oughstun and Sherman [4] to describe the location. Recently, a new 
approximation for this location was obtained in [7]. In this work, however, 
we shall employ a numerical approximation of the saddle point solution ob- 
tained with the help of the Mathematica computer program, and based on 
interpolation techniques. 

As in Oughstun and Sherman study, we deform the original contour of inte- 
gration to the Olver type contour P{9) [8] which passes through the near and 
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distant saddle points. The pole contribution to the asymptotic expansion of 
(9) occurs if in the process of contour deformation one or more poles of 'it^(a;) 
are crossed. From the series representation of the function '^/^(a;) [5] 

U0(u;) = --2il ^ + •••), (10) 

it follows that the integrand in (9) has an infinite set of poles cu — cUc — 2im(5, 
m = 0, 1, 2, • ■ in the half-plane Im u; < 0, which are located along a line, 
parallel to the uo imaginary axis. If j3 is big enough, only the real pole uo = oOc 
is of importance, since the remaining poles are not crossed in the process of 
contour deformation. If, however, /? is small, one or more of the remaining 
poles can be crossed, and their contributions must then be taken into account. 

Let 9s be the value of ^, at which the deformed contour crosses the pole 
&X, uj — ujc in (9), ujc being real and positive. Then, by the Cauchy theorem, 

r 0, e< 0s, 

AJz,t) = { (11) 
I e-f-="^(-=)sin[fo;,(n,(o;,)-^)], 6 > Os- 

Here, nr{uJc) and ni{ujc) stand for real and imaginary parts of n{u!c), respec- 
tively. 

Upon introducing the amplitude attenuation coefficient [4] 

a{ujc) = —UiiuJc), (12) 
c 



and the propagation factor 



C{ujc) = —ririuc), (13) 



c 



Ac{z, t) can be written down as 



0, e<es, 

AJz,t) = { (14) 

g-^a(a;c) gin [((o;,)^ - CU^t], 9 > Os- 

It then follows that for real and positive uJc the pole contribution to the 
asymptotic expansion of A{z, t) oscillates in time at the frequency uJc and 
decreases along its propagation distance z with time independent attenuation 

coefficient a{ujc). 

The pole contribution (14) represents a discontinuous function of 9, while 
the integral representation of A{z, t) changes continuously with 9. As pointed 
in [4] , this fact is of httle significance if z is finite and the pole is bounded away 
from the dominant saddle point at a; = uja- Denote X{uj,9) — Re (f){uj,9). 
Then Q-i^/^)^{^c,d) jg negligible in comparison to the saddle point contribution 
which has the magnitude _ Hence, the discontinuous behaviour 

of Ac{z, t) is then also negligible. 
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4 Uniform asymptotic expression for Ac{z, t) 



The situation becomes different if tlie dominant saddle point approaclies tlie 
pole at a; = ujc- In this case XiujciO) is comparable with X{u}s,0) and so 
are the pole and the branch point contributions to the asymptotic expansion 
of A{z,t). To obtain a continuous asymptotic representation for Ac{z,t), a 
uniform approach, as proposed by Bleistein and Handelsman [10] will here 
be used. 

Let us consider the first pole at a; = uJc- From (9) and (10) it follows that 
the function to be asymptotically evaluated is 

A,{z,t)^-—Re\ du}, (15) 

where \ = z/c. The function (f){uj,6) has a saddle point at c<j = Ug, where 
ujg = ujs{6). (One can verify that the contour P{0) through a near (distant) 
saddle point u = Ug makes the angle 7r/4 (37r/4) with the real axis.) 

In accordance with the Bleistein and Handelsman method we introduce 
a new variable of integration r, defined by 

0(u;,^) = -y-7T + p = *(T,^). (16) 

The quantities 7 i p are chosen so that r = — 7 is the image of the saddle 
point cu — cus and r = is the image of a; = cuc- Then, 

p{9)=(l){uj,,9) and ^(9) = ^2[(l){u;s,e) - (t){u;,,e)]. (17) 

The complex-valued function ^(9) is defined such that it is a smooth function 
of 9 when its argument varies in the interval — tt < Arg 7 < tt. 
One finds from (16) that 



r + ^^^2[<l>(u;s,e)-<l>(u;,e)], (18) 
and hence for cu near cug'. 

T + 7 fti '\/-0a;a;(w8, 0){ljJ - O;,) [1 + 0{uJ - wj]. (19) 

The steepest descent path through the saddle point r = —7 runs parallel to 
the real axis. Upon using (16) in (15) the function Ac{z,t) takes the form 

A,iz,t) = -^Rel[ ^^^e^*(-^)dr|, (20) 
27r [Jc{d) T J 

where 

G„(r,0)=^§^^ (21) 
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and C{9) is the image of P{9) under (16). 
We now expand Gq in the form 

Go(r, 9)=ao + a.r + r(r + 7)i7o(r, 9), (22) 

where Hq{t, 9) is a regular function of r. Since the last term vanishes at both 
critical points r = —7 and r = 0, the coefficients oq and ai are given by 

r(nn^ H go(0,g)-Go(-7,^) 
0'0 = Go{0,9) and Oi = . (23) 

7 

By L'Hospital's rule: 

hm ^ = hm — , (24) 

and hence 

^0(0,^) = !. (25) 

Furthermore, from (19) 

duj 1 
hm — = , 26 

and thus, 

G'o(-7,^) = , \ ■ (27) 

In this manner we obtain 

11 1 

oo = 1 and Oi = — I , . (28) 

If now (22) is inserted into (20), and the resulting canonical integrals ([10]) 
are expressed by special functions, the following result is found 



A,(z,i) = ^Re|e"^ 



W^_i(VA7) + ^WQ{^^\-t) 
V A 



+ i?o(A,^) , (29) 



where, 

Wo(^) = V27re^ and W_i{z) = i H e'^^ds. (30) 
The remainder of the expansion, i?o, is given by 

i?o(A, 9) - A"^ / ^(t, ^)e^*("'^) dr, (31) 

with 

Gi(r,^^) = r^. (32) 
6 



In arriving at (31) we integrated the last term in (22) by parts and neglected 
the boundary contributions as being asymptotically negligible. 

The function W^i can be expressed in terms of the complementary error 
function erfc(z) = 2/{^) e-''' ds. By using (28) and (30) in (29) we 
arrive at the following uniform asymptotic representation 



Ac{z,t) 



Re 



(33) 



2erfc(m'2 




■ CO. 



ujc interacts 



of the main signal in the medium, provided only the first pole a; 
with the saddle point. 

This asymptotic formula applies for any a/A7. In particular, if 7 ^ 0, the 
components of the last parentheses blow up, but their sum remains bounded. 

If yA|7| is large, the error function in (33) can be approximated by its 
asymptotic expansion (comp. [11]) 



eric{iy) = r){y) — 



+ 0{y-') 



(34) 



where 



0, -vr < Arg(y) < 0, 

tliy) = { 1> Arg(y) = -TT or 0, (35) 

2, < Arg(y) < TT. 

Upon using this expansion in (33), the non- uniform asymptotic representa- 
tion of the main signal evolution results: 



le 



A^{z,t) - Re <^ — -77(7) + 



(36) 



If Arg(7) < 0, which occurs when the pole &i uo = oOc is located to the right 
with respect to the contour P{9), the main signal is absent in Ac{z,t), and 
only the term that appears is that proportional to {—(t>ujuj{cOs.,9))~^/'^. This 
term can be interpreted as due to the saddle point ujs- 

If, Arg(7) > 0, which occurs after the contour crosses the pole, then in 
addition to the term proportional to {—(t>ojuj{uJsi ^))~^^^) a new term appears 



'-UJc[nr{uJc) - 0] 



(37) 



It represents the main signal and its form fully agrees with (11). 

In this manner we have obtained both uniform and non-uniform asymp- 
totic representations for the evolution of the main signal in the medium. 
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which are described by (33) and (36), respectively. While the uniform rep- 
resentation applies for any value of \/A7, the non-uniform representation is 
valid only for sufficiently large values of \/A|7|. 

One remark should now be made. The applied theory assumes that the 
saddle point is of the first order, i.e. 0a)a)(a;s, 0) is never zero. In the present 
context this assumption is satisfied everywhere except for the special value 
of = ^^1, where two coalescing near simple saddle points merge on the uj- 
imaginary axis to form a saddle point of the second order. Hence 0^^^ = 
at ^ = ^1, and consequently both asymptotic representations of Ac{z^t), as 
given by (33) and (36), are there invalid. Therefore, strictly speaking, if the 
carrier frequency ujc lies below anomalous dispersion region, (33) is a uniform 
if ^ > Oi. 

5 Numerical example 

A numerical example is now given to illustrate the results obtained in the 
previous section. It is assumed that the Lorentz medium is described by 
Brillouin's choice of medium parameters 

h = X 10^^s-\ LOo = 4.0 X 10^^s-\ S = 0.28 x 10^^s"\ (38) 

and additionally, A = 3.0 x 10"^^, and lOc = 2.0 x lO^^s"^. The latter choice 
implies that in this example the saddle point in question is the near one. 

Let us first suppose that the parameter P in (6) is large enough, say of 
the order of 10^^ or more, to ensure that the second pole lu\,2 = 2^/3 is 
sufficiently distant from the contour P, and, in particular, it is not crossed 
in the process of the contour deformation. Then only the real pole at u = Uc 
is of interest. Under this assumption the real and imaginary parts of 7(6*), as 
given by (16), are shown in Fig. 1. In order to determine numerical values of 
the function u}s{0) an interpolation technique provided by the Mathematica 
computer program has been employed. The evolution of the main signal, 
as predicted by the uniform asymptotic representation (33), is depicted in 
Fig. 2. The anomaly in the plot at 9 = 61 ~ 1.5, results from vanishing 
(f>ujuj{!^sid) ai 9 = 9i, and, as discussed in the previous section, the result 
obtained from (33) breaks down there. Fig. 3 shows the corresponding plot 
obtained from (33), in which the term proportional to l/-\/— has been 
dropped. 

Assume now that P is sufficiently small, such that the second pole at 
i^c2 = <^c ~ 2i/3 can appear close to, or be crossed by the deformed contour P. 
To fix our attention let (3 — 5.0 x lO^^s. In this case an expression similar to 
(33) must be added to the asymptotic representation for Ac{z,t). By virtue 
of (10), the expression should be multiplied by the factor —2, and u!c2 should 
replace lJc- 
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Fig.l i?ea/ and imaginary parts of the function j{0). Here, uJc = 
2.0 X lO^^s"^ and the medium is described by Brillouin's choice of 
parameters. 




Fig. 2 Representation of the main signal in the medium described by 
Brillouin's choice of parameters, based on Eq.(33). Here, uJc — 2.0 x 
W^s-^, A = 3.0x lO-^^s-^ 



9 



Ac ( z , t ) 
0.075 - 




-0.075 



Fig.3 Representation of the main signal evolution in the medium de- 
scribed by Brillouin's choice of parameters, based on Eq.(33) with 
deleted term proportional to {—(j)u)u){<-^s,d))~^^'^- Here, ujc — 2.0 x 
lO^^s-^, A = 3.0x lO-i^s-^ 



Ac (z, t) 

O.lr 




Fig.4 Representation of the main signal in the medium described by 
Brillouin's choice of parameters. Here, uJc — 2.0 x lO^^s"-*^, A = 
3.0 X lO-^^s-i and /3 = 5.0 x 10^1 
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The corresponding plot is shown in Fig. 4. It is seen that now the growth 
of the main signal is slower then in Fig. 2. 

6 Conclusions 

In this paper the problem of electromagnetic signal propagation in a disper- 
sive Lorentz medium is considered. It is assumed that the exciting signal is 
turned on at a finite time instant. The signal rapidly oscillates and its en- 
velope is described by a hyperbolic tangent function. While propagating in 
the medium, the signal splits into three components: Sommerfeld and Bril- 
louin precursors, and the main signal. In this work we find both non-uniform 
and uniform asymptotic representations for the main signal evolution. The 
former representation is readily obtainable by residues. The latter repre- 
sentation is constructed with the help of Bleistein-Handelsman method of 
uniform asymptotic evaluation of integrals with nearby simple saddle point 
and an algebraic singularity. We show, how the uniform representation, ex- 
pressed in terms of complementary error integral, reduces to the non-uniform 
representation. The results here obtained are illustrated with a numerical ex- 
ample. This paper is a complement to our earlier works on Sommerfeld and 
Brillouin precursors ([12], [13]). 
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